Submitted to ApJ on 2008 September 2. Accepted 2009 January 27. 

Preprint typeset using I^T^X style cmulatcapj v. 03/07/07 

DIFFERENT TYPES OF ULTRALUMINOUS X-RAY SOURCES IN NGC 4631 
Roberto Soria 1 , and Kajal K. Ghosh 2 

Submitted to ApJ on 2008 September 2. Accepted 2009 January 27. 

ABSTRACT 

We have re-examined the most luminous X-ray sources in the starburst galaxy NGC 4631, using 
XMM-Newton, Chandra and ROSAT data. The most interesting source is a highly variable supersoft 
ULX. We suggest that its bolometric luminosity ~ a few 10 39 erg s _1 in the high/supersoft state: this 
is an order of magnitude lower than estimated in previous studies, thus reducing the need for extreme 
or exotic scenarios. Moreover, we find that this source was in a non-canonical low/soft (kT ~ 0.1-0.3 
keV) state during the Chandra observation. By comparing the high and low state, we argue that the 
spectral properties may not be consistent with the expected behaviour of an accreting intermediate- 
mass black hole. We suggest that recurrent super-Eddington outbursts with photospheric expansion 
from a massive white dwarf (M wc i > 1.3M@), powered by non-steady nuclear burning, may be a 
viable possibility, in alternative to the previously proposed scenario of a super-Eddington outflow 
from an accreting stellar-mass black hole. The long-term average accretion rate required for nuclear 
burning to power such white-dwarf outbursts in this source and perhaps in other supersoft ULXs is 
« 5-10 x 10~ 6 Af© yr _1 : this is comparable to the thermal-timescale mass transfer rate invoked to 
explain the most luminous hard-spectrum ULXs (powered by black hole accretion). The other four 
most luminous X-ray sources in NGC 4631 (three of which can be classified as ULXs) appear to be 
typical accreting black holes, in four different spectral states: high/soft, convex-spectrum, power-law 
with soft excess, and simple power- law. None of them requires masses > 50Mq. 
Subject headings: X-rays: binaries — X-rays: individual (NGC 4631) — black hole physics 



1. INTRODUCTION 

The most luminous non-nuclear X-ray sources in 
nearby galaxies occur in regions of current or recent star 
formation. Some of them have X-ray luminosities ex- 
ceeding the isotropic Eddintgton luminosity ^Edd f° r an 
« 10M© black hole (BH); they are commonly labelled 
ultraluminous X-ray sources (ULXs). The conservative 
interpretation is that the large majority of ULXs are 
the upper end of the high-mass X-ray binary popula- 
tion, powered by an accreting BH formed from "normal" 
stellar processes. If so, their (apparent) extreme lumi- 
nosity is due to any of the following three reasons, or to 
a combination of them: moderately unisotropic emission 
(King et al. 2001; King 2008); mildly super-Eddington 
luminosities (Begelman 2002, 2006; Ohsuga & Mincshigc 
2007); extremely heavy stellar-mass BHs, with masses 
~ 30-70M Q (Pakull & Mirioni 2002). Those scenarios 
also require mass accretion rates m > 1, where the ac- 
cretion parameter rh = M/A^Edd ~ ^0.1c 2 M^ /L^dd- 

Alternatively, there is still room for the more intrigu- 
ing hypothesis that at least some ULXs are powered by 
intermediate-mass BHs (Miller & Colbert 2004). 

In the absence of direct kinematic measurements (be- 
cause of the faintness of their optical counterparts), X- 
ray spectral and timing studies have been used to try and 
constrain BH masses in ULXs. Such model-dependent 
arguments rely on the (expected) simple scaling of char- 
acteristic variability timescales and disk temperatures 
with BH mass, and on the correspondence of ULX spec- 
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tral states with the "canonical" states of Galactic BHs, 
for which the mass is accurately known. Unfortunately, 
ULXs do not appear to have the same state-transition be- 
havior as Galactic BHs; for example, the most luminous 
sources are rarely found in a high/soft state, dominated 
by a standard accretion disk (Soria & Kuncic 2008). The 
X-ray spectra of some ULXs are dominated by an un- 
broken power law, with photon index T ~ 1.5-2, whose 
physical origin is still unclear. Others have a broad com- 
ponent with a steepening or downward curvature above 
~ 5 keV and sometimes a small soft excess below ^0.5 
keV; this kind of spectrum may come from a slim disk, 
or from the inner region of a standard disk, heavily mod- 
ified by Comptonization when rh > 1. However, there is 
no clear gap between the two kinds of spectra, and the 
phcnomcnological classification of a source in either class 
usually depends on the signal-to-noise ratio available in 
the observations. 

A small subsample of ULXs stands out from this gen- 
eral spectral classification: they have a thermal spectrum 
with temperatures < 0.1 keV and no emission above 1 
keV. This is similar to the spectrum of classical super- 
soft sources in the Milky Way and Magellanic Clouds, but 
their luminosity is one or two orders of magnitude higher. 
The two most luminous supersoft ULXs are M101 ULX-1 
(Kong & Di Stefano 2005) and NGC 4631 XI (Carpano 
et al. 2007). Both sources are strongly variable or tran- 
sient; when in a high state, their blackbody luminosity 
is ~ 10 40 erg s^ 1 . Other supersoft ULXs reaching bolo- 
metric luminosities > 10 39 erg s _1 have been found in 
M81 (Swartz et al. 2002), in the Antennae (Fabbiano et 
al. 2003) and (two) in NGC 300 (Carpano et al. 2006; 
Kong & Di Stefano 2003). Their thermal spectra can, in 
principle, provide tighter constraints on the size of the 
emitting region, and hence more significant tests for the 
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geometry of the accretion flow and the nature and mass 
of the accretor. 

While ULXs with a harder (power-law-like, slim-disk or 
Comptonized) X-ray spectrum may be interpreted as the 
upper end or the natural extension (cither in BH mass or 
accretion rate) of stellar-mass BHs, supersoft ULXs ap- 
pear like the upper end of nuclear-burning white dwarfs, 
which cannot be more massive than w IAM@. Therefore, 
their high apparent luminosities are even more difficult to 
explain. The white dwarf scenario may be salvaged if su- 
persoft ULXs are seen in a transient outburst phase, well 
above the Eddington luminosity of a white dwarf (~ 10 38 
erg s _1 ). The other main competing scenarios for the 
emitting region in supersoft ULXs are a strong outflow 
from a stellar-mass BH accreting at a super-Eddington 
rate (to 3> 1 ), or a standard disk around an intermediate- 
mass BH. For the disk to be so cool (kT in < 0.1 keV), 
the BH mass needs to be > 10 4 M Q . 

In this paper, we re-examine the physical interpreta- 
tion of the supersoft ULX in NGC4631, by comparing 
the XMM-Newton observations taken when the source 
was in a more luminous state, with earlier Chandra ob- 
servations in a lower state. We also discuss the physical 
nature of the other four most luminous X-ray sources in 
the XMM-Newton dataset. 

2. OBSERVATIONS AND DATA ANALYSIS 

NGC 4631 is a late- type starburst galaxy (Hubble type 
SB(s)d), seen nearly edge-on (Figure 1), at a distance of 
7.6 Mpc (Seth, Dalcanton & de Jong 2005). In addi- 
tion to a large number of giant starforming complexes, 
it has one of the best examples of galactic fountains, 
outflows and hot gas above the disk plane (Strickland 
et al. 2004a,b; Wang et al. 2001). From its integrated 
far-infrared luminosity, its star-formation rate is w 3M Q 
yr" 1 (Strickland et al. 2004a; Soifer et al. 1989; Ken- 
nicutt 1998). In the X-ray band, NGC 4631 was stud- 
ied with Einstein (Fabbiano, Kim & Trinchieri 1992), 
ROS 'AT (Vogler & Pietsch 1996; Read, Ponman & Strick- 
land 1997; Liu & Bregman 2005), Chandra (Wang et al. 
2001) and XMM-Newton (Tiillmann et al. 2006a,b; Feng 
& Kaaret 2005; Winter, Mushotzky & Reynolds 2006, 
2007; Carpano et al. 2007). Here, we focus on the five 
brightest point-like sources in the XMM-Newton dataset 
(including four ULXs), and in particular on the supersoft 
ULX (Carpano et al. 2007), whose nature is still contro- 
versial. 

Chandra observations of NGC 4631 with the Advanced 
CCD Imaging Spectrometer (ACIS) were carried out on 
2000 April 16, for 60 ks (ObsID 797; Principal Investiga- 
tor: Daniel Wang). We retrieved the data from the pub- 
lic archive, and analysed them with the locally-developed 
software tool LEXTRCT (Tennant 2006). Source detec- 
tion in LEXTRCT was performed using a circular Gaus- 
sian approximation to the point spread function (PSF), 
which gives higher weight to sources with a central con- 
centration of events. Point-source counts and spectra 
were extracted from within the 95% encircled-energy 
aperture of the model PSF. The background was ex- 
tracted from annular regions surrounding the sources, 
except in crowded regions of the field where we used back- 
ground regions adjacent to the sources. The background- 
subtracted counts within the source regions were scaled 
to obtain the aperture-corrected count values. The 
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Fig. 1. — Multi-band view of NGC 4631. The true-color images 
from top to bottom are: Chandra/ ACIS (red = 0.3—1 keV; green = 
1-2 keV; blue = 2-8 keV); XMM-Newton/EPIC (red = 0.25-1 keV; 
green = 1-2 keV; blue = 2-10 keV), XMM-Newton/OM (red = U 
filter; green = UVW1 filter; blue = UVW2 filter); Spitzer/IRAC 
(red = 5.8 ^m; green = 4.5 /im; blue = 3.6 fim). The five brightest 
sources in XMM-Newton, labelled XI through X5, are the target 
of our study. In all images, North is up and East is left. 
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TABLE 1 

Luminous X-ray sources in NGC 4631 



Source ID ROSAT ID 


R.A. 


Dec. 


Emitted Luminosity ab (erg s x ) 








Chandra/ ACTS XMM-Newton/EPIC ROSAT/RRl ROSAT/PSPC 



XI 

X2 
X3 
X4 
X5 



H13 



H8 
H7 



12 42 15.96 
12 42 11.13 
12 42 06.07 
12 41 57.35 
12 41 55.56 



32 32 49.4 
32 32 35.8 
32 32 46.5 
32 32 03.2 
32 32 16.9 



4.0 



10 37 
x 10 3 



+0.3 v 1n 39 
0.3 



3.9±° 2 2 x 10 38 
3±j x 10 37 



„+0.5 
J -0.9 
1+0.2 

-0.2 

4 0+ 02 

,+0.2 
-0.2 
1+0.2 
-0.2 



3.6: 

3.1 H 



2.r 

5.0" 1 



x 10 39 
x 10 39 
x 10 38 
x 10 39 
x 10 39 



5.7±i;| x 10 39 6.7t\i x 10 39 



9^2 x 10 38 



3.5^ 



x 10 3 



2.9 



+0.4 



x 10 3 



a For the supersoft source XI: bolomctric luminosity (E > 13.5 eV); for the other sources: unabsorbed luminosity in the 0.3-10 keV band. 
b The error ranges listed for Chandra and XMM-Newton luminosities come from our spectral fitting; the error ranges for ROSAT 
luminosities include only the Poisson uncertainty in the HRI and PSPC count rates, after assuming the XMM-Newton best-fitting model 
for counts-to-flux conversion. 



TABLE 2 

Extrapolated bolometric luminosity of XI for different fitting models 



XSPEC model 


xl 


kT hh (cV) 


i-bol (erg 


s- 1 ) 


phabSQ a ] *phabs*bb 


2.20(87.9/40) 




1.5 


X 


10 41 


phabSGal*phabs*(bb+gauss) 


1.40(52.0/37) 




4.5 


X 


1Q 40 


phabSGaj *phabs* (bb+ray ) 


1.80(68.3/38) 




3.6 


X 


10 39 


phabSG a i*phabs* (bb+ray) *zedge 


1.23(44.3/36) 


91+ 7 


3.6 


X 


10 39 


phabSGal*( ra y + ph a bs*bb) 


1.50(57.0/38) 




2.0 


X 


1Q 40 


phabSGal* (ray+phabs*bb) *zedge 


1.24(44.6/36) 


85±% 


5.4 


X 


10 39 



background-subtracted point-source detection limit is 14 
counts for the 2.8 minimum sigal-to-noise ratio (S/N) 
threshold and a minimum 5 a above background. For 
timing analysis, we binned the X-ray light curves of the 
brightest sources into 1000-s bins, and computed \ 2 tests 
against a constant flux hypothesis. For spectral analysis, 
we generated spectral redistribution matrices and ancil- 
lary response files with the Chandra X-ray Center soft- 
ware CIAO version 3.4. We then used XSPEC version 12.0 
(Arnaud 1996) to fit the point-source spectra. 

A 55-ks XMM-Newton observation with the European 
Photon Imaging Camera (EPIC) was carried out on 
2002 June 28 (ObsID 0110900201; Principal Investigator: 
Michael Watson). We downloaded the public- archive 
data and processed them with the XMM-Newton Sci- 
ence Analysis System (SAS) version 6.5.0. We used LEX- 
TRCT for source detection, and standard XMMSELECT 
tasks within the SAS for source and background region 
extraction. The radius of the source extraction regions 
was 20", except for X4, where we used a 15" radius to 
reduce contaminations from the nearby brighter source 
X5. Background extraction regions were chosen around 
the source regions, in a suitable way to avoid contami- 
nation. After building response and ancillary response 
files with rmfgen and arfgen we used XSPEC for spectral 
analysis of the brightest sources (Table 1). To improve 
the signal-to-noise ratio, we co-added the EPIC pn and 
MOS spectra, with suitably averaged response functions, 
using the method of Page et al. (2003). 



To investigate the long-term variability of the ULXs, 
we also re-analyzed the archival ROSAT/RRl and PSPC 
observations carried out between 1991 December and 
1992 December (Voglcr & Pietsch 1996; Read, Ponman 
& Strickland 1997). We applied astrometric corrections 
to the ROSAT data using the Chandra source posi- 
tions. We extracted source counts from circular regions 
of 30"radius, and background counts from source- free cir- 
cular regions of 3'radius. We used WebPIMMS with the 
best-fitting spectral parameters from the XMM-Newton 
study to convert R OS A T source count rates into fluxes. 
Three (X1=H13, X4=H8, and X5=H7; see Table 1 and 
Vogler & Pietsch (1996)) of the five luminous targets of 
this study were also found in ROSAT, with some vari- 
ability over the various exposures. 

3. THE SUPERSOFT ULX 

XI was detected as a luminous supersoft source during 
the ROSAT/PSPC observations of 1991 December 15 - 
1992 January 04 (a total of 18.4 ks; see Fig. 4 and Table 
1 in Vogler & Pietsch (1996)). It was not detected in the 
shorter (3.5 ks) ROSAT/PSPC observations of 1992 May, 
suggesting a count-rate decline by at least a factor of 3. 
It was detected again by ROSAT/RRl in 1992 December; 
however, since the HRI does not provide spectral infor- 
mation, we cannot tell whether it was again in a super- 
soft state. In the Chandra observation from 2000 April 
16, the source was faint and soft. Finally, in the XMM- 
Newton observations of 2002 June 28, XI appeared again 
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TABLE 3 

Best-fit parameters for the coadded EPIC pn and MOS 

SPECTRUM OF XI IN ITS HIGH/SUPERSOFT STATE. SPECTRAL 
MODEL: phabs*phabs* (raymond-smith + bbody) *zedge. VALUES 
IN BRACKETS WERE FIXED. ERRORS ARE 90% CONFIDENCE LEVELS 
FOR 1 INTERESTING PARAMETER (A% 2 = 2.7). 



Parameter XMM-Newton Value 



^H.Gal* 


(1.3 x 10 20 ) 


N n 


2.4+° j x 10 21 


kT rs (kcV) 


44+0-06 
-0.06 


Z(Z Q ) 


(1.0) 


K IB h 


9.2± 2 ; 9 x lO- 6 


kT bb (kcV) 


o.o9i±g:8S 




5.7±° 5 5 x 10-6 


E odge (kcV) 


i m+ ' 01 


T edge 


2 0+ 1 ' 2 



X 2 /dof 1.23(44.3/36) 



/0.3-10 d 


3.3±°1 x 10- 14 


/().3-10,rs C 


x 10 -15 


/i-io f 


1A -ol x !0~ 15 


/l-l(),rs S 


I.OIq J X 10_15 


ibol h 


3.6+° 9 5 x 10 39 


-^bol,rs 


2.6±°j x 10 38 



a From Kalberla et al. (2005). Units of cm- 2 . 

b Raymond-Smith model normalization K IS = 10 -14 /{47r [dj\ (1 + 
z)] 2 } J n e n [fdV, where is the angular size distance to the 
source (cm), n e is the electron density (cm -3 ), and nn is the 
hydrogen density (cm -3 ). 

c Blackbody model normalization -R"bb = L^g/D 2 where L^g is 
the source luminosity in units of 10 39 erg s -1 and Dio is the dis- 
tance in units of 10 kpc. 

d Observed flux in the 0.3—10 kcV band; units of erg cm -2 s -1 . 

°Observed flux in the 0.3—10 keV band, in the Raymond-Smith 
component; units of erg cm -2 s -1 . 

f Observed flux in the 1-10 keV band; units of erg cm -2 s -1 . 

s Observed flux in the 1-10 keV band, in the Raymond-Smith 
component; units of erg cm -2 s -1 . 

h Unabsorbed luminosity for E > 13.6 eV; units of erg s -1 . 

'Unabsorbed luminosity for E > 13.6 eV, in the Raymond-Smith 
component; units of erg s _1 . 

as a luminous supersoft source (Carpano et al. 2007), in 
a similar state to the 1991 detection. Here, we briefly 
summarize the spectral results when the source was in 
a high state, and then compare them with the low-state 
observation. 

3.1. High state 

We re-extracted the XMM-Newton data and coadded 
the pn and MOS spectra with a suitably averaged re- 
sponse function. This is equivalent to fitting them si- 
multaneosly, but provides a better signal-to-noise ratio 
for discrete features. We recover the result of Carpano et 
al. (2007), with a spectrum dominated by a soft black- 



body component, plus residual features (both in emis- 
sion and absorption) especially at m 0.7-1.2 keV (Figure 
2). Such features are already evident in each individual 
EPIC spectrum, as plotted in Figure 4 of Carpano et al. 
(2007), and become more significant when the spectra 
from all three detectors are combined. The visual im- 
pression of such systematic residuals is confirmed by the 
fit statistics: the best-fitting absorbed blackbody model 
has x 2 = 2.20(87.9/40) and can be safely rejected. 

The residual emission and absorption features may not 
seem to affect the bolometric luminosity significantly, 
compared with the dominant blackbody emission. How- 
ever, different models for such residual components have 
the effect of shifting the fitted temperature of the black- 
body component between « 65 and 90 eV, with a dra- 
matic effect on the extrapolated, unabsorbed bolomet- 
ric luminosity. Some examples are summarized in Table 
2. When we model the deviations from a pure black- 
body spectrum with only emission components (e.g., a 
Gaussian as in Carpano et al. (2007), or an optically- 
thin thermal plasma), we find an extrapolated luminos- 
ity w 2-5 x 10 40 erg s _1 . However, we also find that 
there is a statistically-significant absorption feature at 
E = 1.0 ± 0.1 keV. This may be analogous to the Fe-L 
absorption edges found in some Seyfert galaxies (Boiler 
et al. 2003). Ignoring this edge leads to apparently lower 
blackbody temperatures and therefore higher extrapo- 
lated bolometric luminosities. When we include this edge 
in our models (Tables 2 and 3), we obtain bolometric lu- 
minosities as low as « 4-5 x 10 39 erg s _1 , depending on 
whether we assign the same, high intrinsic absorption to 
both the optically-thin and optically-thick thermal com- 
ponents, or only to the latter. On the other hand, the 
isotropic emitted luminosity 2 x 10 39 erg s _1 in the 0.3-10 
keV band provides a solid lower limit to the bolometric 
luminosity. 

In our best-fitting model (Table 3), the optically-thin 
emission component contributes less than 1/10 of the ex- 
trapolated bolometric luminosity, but about 1/4 of the 
observed flux in the full 0.3-10 keV band, about half of 
the unabsorbed luminosity at energies > 0.7 keV, and 2/3 
of the unabsorbed luminosity at energies > 1 keV. This 
explains why modelling such component has a great effect 
on the fit parameters and inferred luminosity. The unab- 
sorbed luminosity of the optically-thin thermal plasma is 
« 2.0 x 10 38 erg s" 1 (bolometric luminosity « 2.6 x 10 38 
erg s _1 ). 

The 1991 ROSAT/PSPC spectrum is also dominated 
by a blackbody component at kT b h J$ 0.1 keV (Read, 
Ponman & Strickland 1997), and a WebPIMMS estimate 
suggests it may have similar luminosity to the 2002 
XMM-Newton spectrum. However, the ROSAT data do 
not have enough counts and spectral resolution to con- 
strain the temperature (and therefore the extrapolated 
emitted luminosity) more accurately. If we assume that 
the spectral model was the same as in the XMM-Newton 
observation (choosing for example the model listed in 
Table 3), and leave only the relative normalization free 
between the two epochs, we find that the two spectra 
are indeed consistent with being very similar (Figure 3), 
with a ROSAT flux normalization 1.33j^Q |q times higher 
than for XMM-Newton. 

The best-fit blackbody temperature (kT bb = 0.09 ± 
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Fig. 2. — Model fits and \ 2 residuals for the combined XMM- 
Newton/EPlC spectrum of XI: simple blackbody (blue curve; 

~ 2.2), and blackbody plus thermal-plasma emission plus ab- 
sorption edge (red curves; xi ~ 1-2; best-fitting parameters in 
Table 3). 
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keV) and 
4 x 10 39 erg s" 1 in 2002, and w 5 x 10 39 
in 1991) correspond to a characteristic radius w 2.1 x 10 9 
cm or w 2.4 x 10 9 cm (in 2002 and 1991, respectively) for 
the optically-thick emitting surface. We should also keep 
in mind that a simple blackbody approximation tends to 
overestimate the true luminosity of supersoft sources by 
a factor of a few (Kahabka & van den Heuvel 1997), so 
the true luminosity may be even lower. 

3.2. Low state 

From the Chandra data in the low state, we found 
w 13.5 net ACIS counts, corresponding to a count rate 
of (2.3 ±0.7) x 10 -4 ct s _1 (a factor of two less than esti- 
mated by Carpano et al. (2007)). The breakdown of the 
net counts in different energy bands is: f» 12.2 counts at 
0.3-1.5 keV, and « 1.3 counts at 1.5-8.0 keV. Dividing 
the interval into three bands, we get: ~ 10 counts at 
0.3-1.0 keV, w 2.5 counts at 1.0-2.0 keV, < 1 counts at 
2.0-8.0 keV. 

So, although there are not enough counts for detailed 
spectral fitting, we have at least a strong indication that 
the source was very soft (or "quasi-soft" in the definition 
of Di Stefano & Kong (2004)). Even if we assume no 
intrinsic absorption, the count distribution rules out a 
power-law spectrum with photon index T < 2 at the 90 
per cent confidence level (using the Cash statistic, Cash 
(1979)). The optically-thin thermal plasma component 
fitted to the high-state spectrum is clearly inconsistent 
with the lower flux detected in the low state (Figure 4). 
More generally, optically-thin thermal plasma models (at 
fixed solar metallicity) are also ruled out at the 90 per 
cent confidence level (the best-fitting model has a Cash- 
statistic parameter = 21.8/14) Soft, optically-thick ther- 
mal emission (blackbody or disk-blackbody) is a much 
better model for the observed count distribution. We es- 
timate a disk-blackbody temperature kT m = 0.22^q' 32 
keV (90 per cent confidence limits; Cash-statistic fit 
parameter — 13.6/14) or a simple blackbody tempera- 
ture fcT bb = 0.18±g;^ keV (Cash-statistic fit parameter 
= 13.2/14). The best-fitting unabsorbed luminosity is 




1 2 
Energy (keV) 

Fig. 3.— Combined fit to the XMM-Nemton/EPIC (solid line) 
and RO/Sj4T/PSPC (dashed line) binned spectra for the supersoft 
source XI in its high state (datapoints and \ 2 residuals). The 
XMM-Newton data are from 2002 June; the ROSAT data from 
1991 December - 1992 January. The main component of the best- 
fitting model (Table 3) is an absorbed blackbody with fcTbb = 
0.08 ± 0.02 keV; the normalization of the ROSAT spectrum is a 
factor Ri 1.3 higher than that of the XMM-Newton spectrum. 

* 10 37 erg s _1 , although the 68 per cent confidence in- 
terval includes values as high as « 10 38 erg s _1 ; but in 
any case, the luminosity is lower than in the supersoft 
state, as intuitively expected from a simple scaling of 
the count rate between the Chandra and XMM-Newton 
observations. From the few detected counts, we cannot 
strongly rule out the alternative possibility that the lower 
count rate in the Chandra observation is also or mostly 
due to a much higher intrinsic absorption (iVn > 10 22 
cm~ 2 ) of the same, luminous super-soft component. But 
in the absence of independent evidence for that dramatic 
change in absorption, and by analogy with the behaviour 
of other accreting systems (including the supersoft source 
in M101 mentioned earlier, Kong & Di Stefano (2005)), 
here we consider the low-state scenario as the most plau- 
sible. 

The absence of the optically-thin thermal plasma com- 
ponent in low state is somewhat puzzling. The rapid 
change between 2000 and 2002 suggests that it was not 
due to diffuse hot gas at large distances from the source, 
but was instead directly associated with the high state 
or outburst. The PSF and source extraction region are 
of course much larger for XMM-Newton, which may sug- 
gest a contamination from Galactic-scale diffuse emis- 
sion; however, we tested this possibility using the Chan- 
dra images, and we do not find evidence of local enhance- 
ments in the diffuse soft emission. Different choices of 
background extraction regions in XMM-Newton do not 
remove this component. In conclusion, we suggest that 
the most likely explanation at this stage is that the line 
emission and absorption edge are really associated with 
the compact source. If the emission is due to an ex- 
panding BH wind or white dwarf photosphere (scenarios 
outlined in Sections 4.2 and 4.3) we speculate that this 
may be evidence of an optically thick and optically thin 
component in the outflow. 

4. PHYSICAL INTERPRETATION OF XI 

Three scenarios have been considered for this source 
(and for the handful of similar ones in other nearby galax- 
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Fig. 4. — Datapoints from the Chandra observation in the low 
state (each bin corresponding to a single detected count) com- 
pared with the expected contribution from the optically-thin ther- 
mal component alone, using the best-fitting parameters of the high 
state (Table 3). 

ies), based on the properties of its high/supersoft state 
(Carpano et al. 2007): standard disk emission from an 
intermediate-mass BH; super-Eddington outflows from a 
stellar-mass BH; nuclear-burning white dwarf. We dis- 
cuss here how the low-state appearance provides stronger 
constraints. 

4.1. Intermediate-mass BH 

In the standard disk model, the BH mass can be ex- 
pressed as a function of inner-disk temperature and bolo- 
metric disk luminosity (or diskbb normalization constant): 

M ~ Tr^L 1 ^ (Makishima et al. 2000). Assuming 
i?in = 6M and imposing Ah > 1.3 x 10 20 cm~ 2 (Galactic 
line of sight), we used the Cash-statistics fit to the Chan- 
dra data, and we derived a contour plot of the acceptable 
region in the BH-mass versus disk-luminosity space (Fig- 
ure 5). We do not have enough net counts to determine 
the 90% confidence contours, but we can at least iden- 
tify a "most plausible" (68% confidence contours) region 
in that plane. The large uncertainty in the inner-disk 
temperature (from ~ 0.1 to ~ 0.5 keV) is reflected in a 
large range of masses and unabsorbed luminosities. Re- 
gardless of BH mass, the luminosity appears to be al- 
ways < O.OliEdd- However, there are observational and 
theoretical arguments against a disk-dominated low/soft 
state at L < O.OlLEdd- At those luminosities, accret- 
ing BHs are generally found in a power-law-dominatcd 
low/hard state, which is ruled out in this case. Hence, we 
suggest that the low/soft appearance of the source in the 
Chandra observation does not favour the disk-blackbody 
model. 

Moreover, the fitted inner-disk temperature in the 
low/soft state appears to be similar or slightly higher 
than the thermal temperature in the high/supersoft state 
(Carpano et al. (2007), and Table 2). If the thermal 
emission is due to an accretion disk, we expect the disk 
to be cooler, when the net count rate (a proxy for the 
accretion rate and emitted luminosity) is two or three 

orders of magnitude lower; T- m ~ £^{ s 4 k in a standard 
disk. A transition between high/supersoft and low/soft 



states is inconsistent with the well-studied behaviour of 
BH accretion disks. In conclusion, the comparison of 
low- and high-state spectral data does not favour the 
intermediate-mass BH scenario. 

4.2. Stellar-mass BH 

The high/supersoft state is consistent with ther- 
mal emission scattered and collimated by a massive, 
optically-thick disk outflow, launched at the spheriza- 
tion radius around a stellar-mass BH, when the accretion 
rate exceeds the Eddington limit (Shakura & Sunyaev 
1973; Poutanen et al. 2007; King 2008). We have already 
noted that in the Chandra observation, the emitted lu- 
minosity may be two or three orders of magnitude lower. 
Hence, we do not expect the super-Eddington outflow to 
be present in that lower state. 

In the canonical scheme of BH accretion states, inferred 
mostly from the study of stellar-mass Galactic BH bina- 
ries, moderately active (sub-Eddington) BHs arc cither 
in the high/soft state (dominated by a disk-blackbody 
component with fcT; n « 1 keV), or in the low/hard 
state (power-law component with V ~ 1.5-2). Neither 
state is consistent with the observed Chandra spectrum. 
And conversely, no stellar-mass BHs have been observed 
in a disk-dominated state with kT m ~ 0.2 keV and 
Lx ~ 10 37 -10 38 erg s~ x . If the supersoft component 
in the higher state is attributed to a disk outflow during 
super-critical accretion, we would have to conclude that 
the lower-state spectrum does not look like a canonical 
state for stellar-mass BHs. 

An alternative possibility we should consider is that 
the sequence of high and low states is not due to accre- 
tion state transitions, but to a long-term precession of 
the binary system, such that the moderatlcy collimated 
outflow moves in and out of our line of sight. When we 
are looking down the outflow, we may be observing the 
undisturbed, cooler (soft spectrum) outer disk, at radii 
larger than the spherization radius. However, a diffi- 
culty of this interpretation is that the predicted isotropic 
luminosity of the outer standard disk down to the spher- 
ization radius should be « i^Edd ~ 10 39 erg s _1 for a 
stellar-mass BH. This is much higher than observed. 

In conclusion, we suggest that the stellar-mass BH sce- 
nario, although still viable, has not yet provided a per- 
fectly self-consistent interpretation for this source, or at 
least requires a new kind of accretion-state behaviour 
or accretion-disk structure, so far unobserved in stellar- 
mass BHs. 

4.3. Nuclear-burning white dwarf 

The observed temperature and luminosity of XI in the 
low/soft state is consistent with the thermal emission 
from surface hydrogen burning on a massive white dwarf, 
Mwd ~ 1.3-1.35 (Figure 6). Such process can occur at 
accretion rates just below the steady burning rate (~ 
a few 1O" 7 M yr" 1 : Starrfield et al. (2004); Kahabka 
(2004) . The white dwarf is sufficiently hot that hydrogen 
burns on its surface immediately as it is accreted; for this 
reason, the system does not go through classical-nova 
outbursts. A blackbody luminosity « 10 38 erg s _1 at 
a temperature fclbb ~ 0.1 keV are well within the 68% 
confidence limit of the Chandra detection. These values 
correspond to a radius ~ 3000 km for a spherical emitter, 
consistent with the radius of a 1.35-M Q white dwarf. 
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Fig. 5. — Contour plot for the 68% confidence region in the 
BH mass versus bolometric disk luminosity plane, for XI in its 
low/soft state (Chandra observation), if its spectrum is fitted with 
an absorbed disk-blackbody. The dash-dotted line corresponds to 
^disk = O-Ol^Eddi a conventional threshold below which the disk 
is not expected to be the dominant emitting component in ac- 
creting BHs (and therefore a disk-blackbody model would not be 
self-consistent). 
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Fig. 6. — Contour plot for the 68% confidence region in the 
BH mass versus bolometric luminosity plane, for XI in its low/soft 
state (Chandra observation), if its spectrum is fitted with an ab- 
sorbed blackbody. The dashed lines correspond to the predicted 
effective temperatures and luminosities for massive white dwarfs 
(1.25M0 and 1.35Mq) during phases of surface hydrogen burn- 
ing, for different values of mass accretion rates, increasing from 
1.6 X 1O- 8 M yr" 1 to 8.0 X 10- 7 M Q yr" 1 (see Starrfield et al. 
(2004) for detailed explanations of those sequences). 

However, the high/supersoft phases cannot be ex- 
plained with steady surface burning. It was found from 
numerical simulations (Starrfield et al. 2004) that steady 
surface hydrogen burning on a 1.35-Mq white dwarf can 
occur only for accretion rates < 10 _6 M Q yr^ 1 ; above 
that limit, the photosphere of the white dwarf expands 
to large radii and shuts off accretion. This may be what 
is happening in XI, with a sequence of outbursts in be- 
tween phases of more steady surface burning (at reduced 
accretion rate) or simply of surface cooling. The esti- 
mated radii (a factor of 10 higher than in the low/soft 
state, Section 3.1) and blackbody temperatures (a factor 



of two lower) fitted to the high/supersoft spectrum are 
consistent with photospheric expansion. 

From the sequence of X-ray observations between 1991 
and 2002, we speculate that XI is in a transient supersoft 
state about half of the time. Since the emitted luminos- 
ity in the high/supersoft state is ~ a few 10 39 erg s" 1 , 
and that in the low/soft state is ~ a few 10 37 erg s _1 , the 
long-term average luminosity may be ~ 2 x 10 39 . This 
clearly raises two problems, related to the average long- 
term accretion rate and to the peak luminosity. Since all 
the power must ultimately come from nuclear burning, 
the long-term average luminosity requires an average ac- 
cretion rate ~ (0.5-1) x 10~ 5 Afo yr . From stellar evo- 
lution models (Rappaport, Podsiadlowski & Pfahl 2005), 
it appears that even late-type B stars (likely progeni- 
tor of the putative massive white dwarf and likely donor 
star in the system) can provide such extreme rates dur- 
ing thermal-timescale episodes of mass transfer, lasting 
~ 10 5 yr. If the same accretion rate was used to power an 
accreting BH, it would produce a bolometric luminosity 
~ (3-6) x 10 40 erg s _1 , which is similar to the maxi- 
mum luminosity inferred for ULXs; the X-ray luminosity 
distribution of ULXs has a downturn at Lx ~ 2 x 10 40 
erg s _1 (Grimm, Gilfanov & Sunyaev 2003; Swartz et 
al. 2004). In other words, explaining the most luminous 
supersoft sources ever found as nuclear-burning sources 
with average long-term luminosity ~ 10 39 erg s _1 re- 
quires mass-transfer conditions similar to those needed 
to explain the most luminous ULXs as accreting BHs 
with luminosities ~ a few 10 40 erg s _1 . We speculate 
that this may be more than a coincidence. Moreover, 
the alternative ULX scenario of a super-Eddington out- 
flow from a stellar-mass BH (Section 4.2; Carpano et al. 
(2007); King (2008)) also requires mass transfer rates of 
up to » 10 -5 M© yr -1 on the thermal timescale, possibly 
from a B-type donor star. Therefore, the required limits 
on the long-term average mass transfer rate are similar 
for both models. 

As for the peak luminosity ~ a few 10 39 erg s" 1 , this 
is highly super-Eddington for a white dwarf, and a fac- 
tor of 10 above the luminosity produced by steady sur- 
face burning. However, steady burning cannot persist at 
accretion rates w 1O _5 M0 yr -1 . We speculate that the 
high/supersoft phases may be transient super-Eddington 
events (fireball scenario), after which the photosphere 
shrinks again to the white-dwarf surface, and accretion 
resumes. It was estimated (Starrfield et al. 2004) that 
the hydrogen layer involved in surface nuclear burning 
has a mass ~ 10 -6 M Q . Therefore, even if further hy- 
drogen accretion is shut off (i.e., during the transient 
super-Eddington outburst), simply the complete burn- 
ing of this layer can provide a luminosity > 10 39 erg s _1 
for several weeks, or > 10 40 erg s" 1 for several days, 
and the layer itself can be replenished in ~ 1 month, 
before another outburst. A similar process may be driv- 
ing the outburst cycles of the supersoft ULX in M101 
(Kong & Di Stefano 2005), which shows varying tem- 
peratures between sa 50 and 150 eV corresponding 
perhaps to phases of photospheric expansion and con- 
traction. Super-Eddington outbursts powered by non- 
steady episodes of nuclear burning have been observed 
in some Novae, most notably LMC 91, which peaked at 
w 2.6 x 10 39 erg s" 1 (Schwarz et al. 2001). 

In conclusion, we suggest that a fireball white-dwarf 
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model is still a viable scenario for this extreme source 
(and perhaps also for the whole class of supersoft ULXs), 
considering that its true bolomctric luminosity is likely to 
be an order of magnitude less than originally estimated 
(Carpano et al. 2007), and that the required accretion 
rate ~ 10 _5 M Q yr _1 is similar to the rate invoked for 
super-Eddington outflows in the most luminous ULXs. 

4.4. Time variability 

From a timing study of the XMM-Newton observa- 
tion, it was found (Carpano et al. 2007) that the X-ray 
emission has a modulation with an apparent period of 
« 4.1 hr. It is not clear what can produce this period or 
timescale. If it was the binary period, from the period- 
density relation (Warner 1995) we infer an average den- 
sity w 6 g cm -3 inside the Roche lobe of the donor star 
(the mean solar density is « 1.4 g cm -3 ). This is con- 
sistent with a main-sequence star with a mass « 0.5M©. 
Other well-known supersoft sources such as CAL87 are 
known to have a low-mass donor with a binary period 
of a few hours (Callanan et al. 1989). However, for 
NGC 4631 XI, such a low-mass donor is inconsistent with 
the required mass transfer rates in both the stellar-mass 
BH and white dwarf scenarios, and is also at odds with 
the young age of the stellar population in this starburst 
galaxy (i.e., the galaxy is more likely to contain bright 
high- mass X-ray binaries). A B-type donor star is con- 
sistent with a massive white-dwarf compact object, both 
from theoretical arguments (Iben & Tutukov 1994) and 
observationally (Bcrghofer et al. 2000); the lifetime of the 
B-type progenitor of a 1.3M Q white dwarf is w 50 — 70 
Myr. 

If the 4-hr X-ray modulation is not the binary period, 
then it could be due either to a disk precession in the BH 
scenario, or to the rotational period of the white dwarf, 
or to pulsations in the donor star which affect the rate of 
mass transfer. A study of these scenarios is beyond the 
scope of this work. We just point out that 4 hrs is the 
characteristic pulsation period of the B-type /3-Cephei 
stars (Stankov & Handler 2005). Pulsations of a /3-Cep 
donor have been invoked in the past as a possible cause of 
X-ray periodicities in some accreting binaries (Bcrghofer 
et al. 2000; Finley et al. 1992). 

5. OTHER LUMINOUS X-RAY SOURCES 

5.1. X2 

Apparently coincident with a young star cluster, this 
highly-absorbed ULX (Ah ~ 2 x 10 22 cm -2 ), was not 
detected in any ROSAT observation (not surprisingly). 
This source can be classified as a "convex-spectrum" 
ULX (using the terminology of Makishima (2007)). The 
spectral curvature can be formally modelled with a stan- 
dard disk-blackbody spectrum (Figure 7 and Table 4), 
with a color temperature kT m w 1.3-1.5 keV (in the 
XMM-Newton and Chandra data, respectively). How- 
ever, as is generally the case in this class of ULXs, such 
temperatures are too high for the estimated luminosities 
(« 3 and 4 x 10 39 erg s _1 , respectively). One possibility 
is that the emission comes from a slim disk (Watarai, 
Mizuno & Mineshige 2001) rather than a standard disk; 
if so, the mass accretion rate may be an order of mag- 
nitude above the Eddington rate, while the luminosity 
may be ~ LEdd, and the BH mass ~ 20-30M© (for com- 
parison with other sources in a similar state, see Fig. 3 in 



TABLE 4 

Best-fitting parameters for the Chandra/ ACIS and 
XMM- Newton/ EPIC spectra of X2. Spectral model: 
wabs*wabs*diskbb. Values in brackets were fixed. Errors 

ARE 90% CONFIDENCE LEVELS FOR 1 INTERESTING PARAMETER 

(A X 2 = 2.7). 



Parameter 



Chandra Value XMM-Newton Value 



^H,Gal a 


(1.3 x 10 2() ) 


(1.3 x 10 2() ) 


N H 


28.3+J'j x 10 21 


26.4+J'j x 10 21 


kT dhh (keV) 


1 49+0-22 


1 26+ ' 11 


i^dbb b 


5.8+Jj x !0 -3 


8 - 5 -t'.8 X 10_3 



X 2 /dof 


0.65(42.4/65) 


0.78(42.1/54) 


/0.3-10 C 


3.010 2 x 10~ 13 
4.0±°T 3 x III-" 


2.1±°-\ x 10 -13 
3.1±°,; 2 x 10 s9 



a From Kalberla et al. (2005). Units of cm 2 . 

b ^dbb = [Rin(krn)/c((10 kpc)] 2 X cos# where R ln is the apparent 
inner-disk radius and 8 the viewing angle; 8 = is face-on. 
c Observed flux in the 0.3-10 keV band; units of erg cm~ 2 s — 1 . 
d Unabsorbcd luminosity in the 0.3-10 keV band; units of erg s -1 . 
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Fig. 7. — XMM-Newton/EPIC spectrum of X2 (datapoints and 
X 2 residuals), fitted with a highly absorbed disk-blackbody model; 
see Table 4 for the best-fitting parameters. 

Makishima (2007)). Alternatively, the convex spectrum 
may be modelled equally well with a Comptonized com- 
ponent, arising from a warm (kT ~ 3-4 keV) corona (see 
a comparison between the two scenarios in Stobbart et 
al. (2006)). Other purely phenomenological models such 
as a broken power-law (breaking at ~ 4 keV) also provide 
good fits. 

5.2. X3 

This X-ray source is also well modelled with a disk- 
blackbody spectrum (Figure 8 and Table 5), with fcTin « 
1.2-1.4 keV (in the XMM-Newton and Chandra data, 
respectively), similar to the parameters found for X2. 
However, its emitted luminosity in the 0.3-10 keV band 
is only w 4 x 10 38 erg s _1 , constant between Chandra and 
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Fig. 8. — XMM-Newton/EPIC spectrum of X3 (datapoints and 
X 2 residuals), fitted with a disk-blackbody model; see Table 5 for 
the best-fitting parameters. 

XMM-Newton and below the ROSAT detection limit. 
Therefore, X3 is consistent with a stellar-mass BH near 
the upper- luminosity end of its high/soft state. Assum- 
ing a spectral hardening factor « 1.7-2 (Shimura & Taka- 
hara 1995; Gierlihski & Done 2004; Shafee et al. 2006), 
the standard-disk temperature-luminosity relation (Mak- 
ishima et al. 2000) suggests a BH mass w 5-7M©; the 
fitted (apparent) inner-disk radius R ln ~ 30(cos#) -0 ' 5 
km, as expected. 

5.3. X4 

This ULX shows a state transition between the Chan- 
dra and XMM-Newton observations. In Chandra, we 
detect a faint, soft source (0.3-8 keV count rate = 
(1.44±0.16) x 10~ 3 counts s" 1 ), well fitted (xl = 0.88) by 
an optically-thin thermal-plasma model (Table 6) with a 

temperature kT rs = 1.2+Q2 keV. Every other spectral 
model (power-law, disk-blackbody, or any Comptoniza- 
tion models) yields xt> ^ 1-5- The emitted luminosity for 
the Raymond-Smith thermal-plasma model is w 3 x 10 37 
erg s _1 . In the XMM-Newton observations, there is a 
stronger, harder source at the same position, with a 
broad spectrum well fitted by a power-law plus disk- 
blackbody model (Figure 9 and Table 6). The power- 
law index is T w 1.9; there is no evidence of a steepen- 
ing break or spectral curvature near or above ~ 5 keV, 
unlike what we noted for X2. The disk-blackbody com- 
ponent has a color temperature kT m = 0.20+o 05 keV. 
The unabsorbed luminosity in the 0.3-10 keV band is 
(2.1 ± 0.2) x 10 39 erg s _1 , « 30 per cent of which in 
the disk-blackbody component. It is possible that the 
thermal-plasma component seen in the Chandra dataset 
is also present in the XMM-Newton spectrum, although 
we cannot place strong constraints on it (Table 6). 

In summary, X4 is a transient "power-law" ULX (us- 
ing again the spectral classification of Makishima (2007)) 
with a thermal soft-excess at low energies. A physical in- 
terpretation for this class of ULXs is that the standard 
optically-thick accretion disk is directly visible at large 
radii but is replaced or covered by a scattering-dominated 
region (producing a broader, power-law-like spectrum) 
at small radii (see, e.g., the review by Soria & Kuncic 
(2008) and references therein). The apparent inner-disk 



TABLE 5 

Best-fitting parameters for the Chandra/ ACIS and 
XMM-Newton/ 'EPIC spectra of X3. Spectral model: 
wabs*wabs*diskbb. Values in brackets were fixed. Errors 

ARE 90% CONFIDENCE LEVELS FOR 1 INTERESTING PARAMETER 

(A X 2 = 2.7). 



Parameter Chandra Value XMM-Newton Value 



W H ,Gai a 


(1.3 x 10 2() ) 


(1.3 x 10 2() ) 


N H 


2.0±J- °, x 10 21 


i.i±S:| x io 21 


kT dhh (keV) 


I oy + 0.32 

1 -°'-0.21 


1 20+ ' 17 


i^dbb b 


7.8+lj x 10- 4 


13. 8^4 4 x 10~ 4 


x 2 /dof 


0.81(17.9/22) 


0.75(30.7/41) 



/o.3-io c 4.7l ( ,!;i x 10- 14 5.0±°,;i x 10- 14 
Lo.3-io d 3.9«; 2 x 10 38 4.0«; 2 x 10 38 



a From Kalberla et al. (2005). Units of cm" 2 . 

b ^Qbb = [-Rin(km)/(2(10 kpc)] 2 X cos# where R ln is the apparent 
inner-disk radius and 8 the viewing angle; 8 = is face-on. 
c Observed flux in the 0.3-10 keV band; units of erg cm~ 2 s _1 . 
d Unabsorbcd luminosity in the 0.3—10 keV band; units of erg s — 1 . 

radius (which we may identify with the transition ra- 
dius between standard disk and Comptonizing regions) 
i?i n « 1600(cosf5)~ - 5 km. From the relative contribu- 
tion of thermal and (less radiatively efficient) power-law 
components, we speculate that this transition radius is 
> a few times the innermost stable orbit; hence, the BH 
mass is likely to be < 50M Q . 

Evidence of variability on monthly timescales for X4 
was already found in the series of R OS A T observations 
(Vogler & Pietsch 1996). We searched for short-term 
variability in the XMM-Newton observation (Figure 10), 
and found that a constant count rate is statistically ruled 
out (xl = 117.2/83). However, we found no dominant 
frequency or other spectral features in its power-density 
spectrum. 

5.4. X5 

This ULX was in a luminous state in five of the six 
R OS A T observations (Vogler & Pietsch 1996). We found 
it again in a bright state during the Chandra and XMM- 
Newton observations (Figure 11 and Table 7). In both 
datasets, the X-ray spectra are well fitted by a sim- 
ple power-law of photon index r ~ 2. There is no 
statistically-significant evidence of either a soft thermal 
component or a high-energy break or downward curva- 
ture in the power-law. The unabsorbed isotropic lumi- 
nosity in the 0.3-10.0 keV band is (3.8 ± 0.1) x 10 39 
ergs s^ 1 {Chandra) and (5.0 ± 0.2) x 10 39 ergs s" 1 
XMM-Newton, only slightly higher than the range of esti- 
mated luminosities during the ROSAT observations (Liu 
& Bregman 2005), when the same Chandra or XMM- 
Newton model is applied to the ROSAT data 3 . We did 
not find any variability within the individual Chandra 

3 The ROSAT/THil luminosity estimated by Liu & Bregman 
(2005) is a factor of 2 lower than our estimate because they assumed 
only line-of-sight absorption. 
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TABLE 6 

Best-fit parameters for the Chandra/ ACIS and 
XMM-Newton/EPIC spectra of X4. Spectral model: 
wabs* [(wabsi *raymond-smith) + wabS2*(power- law+diskbb)] . 

Values in brackets were fixed. Errors are 90% confidence 
levels for 1 interesting parameter (Ax 2 = 2.7). 



Parameter 



Chandra Value XMM-Newton Value 



JV H,Gal 


^±.0 X _1_U J 




iV H ,i 


3.2li° 2 2 x 10 20 


(3.2 X 10 2() ) 






3.6±°;g x 10 21 


fcT rs (keV) 


1 24+ 053 


(1.24) 


Z(Z ) 


(1.0) 


(1.0) 


iv rs b 


2.7±\ l x 10~ 6 


2.2±£_l x 10~ 6 


r c 




1 88+ ' 13 
i -°°-0.12 


iV P i d 




3.5±g;| x 10 -5 


fcT dbb (keV) 




20+ 11 
u - zu -0.05 


^dbb c 




A k+2.7 

4.5_ 44 


X 2 /dof 


0.88(4.4/5) 


0.98(67.7/69) 



/o.3-10 
io.3-10 g 



O.03±2;{£ x 10 39 



-1.0 

2.1±° 2 x 10 39 



a From Kalberla et al. (2005). Units of cm- 2 . 
b Raymond-Smith model normalization K IS = 



10 -14 /{4ttKi(1 + 
z)] 2 } J n e nifdV, where d\ is the angular size distance to the 
source (cm), n e is the electron density (cm -3 ), and nu is the 
hydrogen density (cm -3 ). 
c Photon index. 

d Units of photons kcV~ 1 cm -2 s — 1 , at 1 keV. 

c -^"dbb = [-Ri n (km)/d(10 kpc)] 2 x cos 9 where Rj n is the apparent 
inner-disk radius and 9 the viewing angle; 9 = is face-on. 
'Observed flux in the 0.3—10 keV band; units of erg cm~ 2 s~ 1 . 
g Unabsorbcd luminosity in the 0.3—10 keV band; units of erg s^ 1 . 

and XMM-Newton exposures. In summary, X5 appears 
to be a typical power-law ULX, perhaps powered by ac- 
cretion onto a BH with a mass < 50M Q . 

6. CONCLUSIONS 

We have studied the nature of the five brightest sources 
in NGC4631, using XMM-Newton and Chandra data. 
Four of them can be classified as ULXs. The most pe- 
culiar ULX, which we label XI, was previously studied 
by Carpano et al. (2007) and identified as a variable su- 
persoft source with an apparent bolometric luminosity 
~ 3 x 10 40 erg s _1 . We re-examined the spectral data 
and found that in fact, its most likely luminosity may be 
only w 4 x 10 39 erg s _1 ; in fact, it could be even lower, 
if we consider that a blackbody approximation tends to 
overestimate the luminosity of supersoft sources. This 
reduces the need for an intermediate-mass BH or other 
exotic scenarios. 

We found that when the source is in a low state ( Chan- 
dra observation), it appears soft (but not supersoft), con- 
sistent with a thermal spectrum at a temperature ~ 0.1- 
0.3 keV, and a luminosity 10 37 < L bo i < 10 38 erg s" 1 . 
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Fig. 9. — XMM-Newton/EPIC spectrum of X4 (datapoints and 
X 2 residuals), fitted with a disk- blackbody plus power- law model; 
see Table 6 for the best-fitting parameters. 
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Fig. 10. — XMM-Newton/EPIC lightcurve of X4, and corre- 
sponding background count rate, showing significant (aperiodic) 
short-term variability. 
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Fig. 11. — XMM-Newton/EPIC spectrum of X5 (datapoints and 
X 2 residuals), fitted with a simple power-law model; see Table 7 
for the best-fitting parameters. 
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TABLE 7 

Best-fitting parameters for the Chandra/ ACIS and 
XMM-Newton/EPIC spectra of X5. Spectral model: 
wabs*wabs*power-law. Values in brackets were fixed. 
Errors are 90% confidence levels for 1 interesting 
parameter (Ax 2 = 2.7). 



Parameter 


Chandra Value 


XMM-Newton Value 


^H.Gal" 

r b 


(1.3 x 10 20 ) 

2.0t'o 2 x 1021 

1 7g +0.08 
^-0.08 

8.6to_7 x 10~ 5 


(1.3 x 10 2() ) 
2.5I°; 2 x 10 21 

9 i 9 +0.06 
z - ±z -0.03 

13.7±"[ 8 7 x lO" 5 


X 2 /dof 


0.95(140.0/147) 


1.13(187.1/166) 


/o.3-l() d 
£o.3-10 C 


4.2+ ( ^;| x 10- 13 
3.8l°;i x 10 39 


4.4+ ( ^J x 10- 13 
5.0±°; 2 x 10 39 



a From Kalberla et al. (2005). Units of cm- 2 . 
b Photon index. 

c Units of photons keV - 1 cm~ 2 s — 1 , at 1 keV. 

d Observed flux in the 0.3—10 keV band; units of erg cm~ 2 s~ 1 . 

°Unabsorbed luminosity in the 0.3-10 keV band; units of erg s -1 . 

We argued that this is inconsistent with an intermediate- 
mass BH. It is also unsual for a stellar-mass BH, which at 
those luminosities is expected to be in a power-law dom- 
inated low/hard state, or in a disk-dominated high/soft 
state with T; n ~ 0.5-1 keV, based on our current knowl- 
edge of canonical accretion states. As an alternative, we 
suggest that transient supcr-Eddington outbursts (fire- 
ball or photospheric expansion) powered by non-steady 
nuclear burning on the surface of a massive white dwarf 
could be a viable scenario, as an extreme subclass of 
supersoft sources. Outbursts due to photospheric ex- 
pansion are expected when the accretion rate exceeds 
w 10~ 6 M Q yr _1 . Based on the sequence of available 
observations from 1991-2002, the long-term average lu- 
minosity of the system is « (1-2) x 10 39 erg s _1 , which 
requires average accretion rates ~ (5-10) x 10~ 6 M Q yr _1 . 
Although very high, such rates are achievable during 
phases of thermal-timescale mass transfer in B stars, and 
are similar to the rates required to explain the most lu- 
minous ULXs powered by BH accretion. Hence, we spec- 
ulate that transient outbursts in nuclear-burning, mas- 
sive white dwarfs may also explain the few other super- 
soft ULXs (all highly variable) found in nearby galaxies. 
Some of those sources would be seen as quasi-soft sources 



(in the definition of Di Stefano & Kong 2004) when they 
are in a low state. 

The origin of the 4-hr X-ray variability remains un- 
explained, whatever the nature of the compact object. 
Given the high mass transfer rate, we would expect a B- 
type donor star filling its Roche lobe; however, a 4-hr pe- 
riod does not allow for massive donors. Alternatively, the 
variability may be due to an accretion disk precession, or 
the rotation of the white dwarf, or /3-Cephei pulsations 
in the donor star. 

The other four brightest sources in NGC4631 are 
almost certainly bone fide BHs, in different accretion 
states. X2 (L x « 3 x 10 39 erg s^ 1 in the XMM-Newton 
observation) is a highly absorbed "convex-spectrum" 
ULX; its X-ray spectrum may be interpreted as emission 
from a slim disk, or from a low-temperature (a few keV) 
Comptonizing region. X3 (Lx ~ 4 x 10 38 erg s _1 ) is a 
stellar-mass BH in its classical disk-dominated high/soft 
state. X4 (Lx ~ 3 x 10 39 erg s -1 ) is a transient power- 
law ULX with a soft-excess at kT rts 0.20 keV; for this 
class of objects, we may be seeing a standard disk outside 
a transition radius, completely replaced or covered by a 
Comptonizing region at smaller radii. X5 (Lx ~ 5 x 10 39 
erg s _1 ) is a pure power-law ULX, with no evidence for 
low-energy soft excess or high-energy steepening. All 
of them are consistent with accreting stellar-mass BHs, 
with masses < 50M Q . However, the relation between 
the different phenomcnological states (high/soft, convex- 
spectrum, power-law, power-law with soft excess) , and in 
particular whether those states are uniquely a function of 
the normalized accretion rate, remains a topic for further 
theoretical and observational investigations. 
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